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▲  SIMPLIFIED  STUDY  OF  OHN-CYCLB  MAfflffiTOHHStODYKAMIC  POWER 
GENERATION  AND  AN  EXAMINATION  OF  THE  POSSIBLE  ROIE 
OF  SUPERSONIC  COMBUSTION  IN  THE  CXCUf . 


|5\  J.  G.  Woodley 


'This  Report  present 


.l-.'-U.  1  ■  ' 


^This^epart  preaenta  a  preliminary  simplified  study  of  sane  of  the 
facets  Involved  in  obtaining  Open-cyal#  M.H.D.  power  generation*  Various 
parts  of  an  idealised  cycle  are  studied  in  detail,  and  examination  has  been 
made  of  the  possibility  of  supersonic  combustion  farming  a  natural  foatiye  of 
M*H*D*  power  generation  cycles  in  which  a  supersonic  M.H.D. .  gener&tingjfauct  is 
used  .  % 


L  In  Appendices  A  to  £,  far  completeness,  formulas  relevant  to  the  various 
stages  of  the  cycle  are  listed  or  derived*. 
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1  INTRODUCTION  -  IDEALISED  M.H.D.  PDfBR  GENERATION  CYC1B 

Fig.  1(a)  shows  an  idealised  M.H.D.  power  generation  cyole,  in  which  the 

_  e  .  .  a  a  j  j»a..  .  ~ n  n.  * ««.*  m  m  m  im  f Ar^tWa  4* sa 

compressor,  nosale  ana  axuusei'  oi  uia  a  jo  w«  ou.o  osomuw* 
is  on  tropic  ally.  In  this  cycle  it  is  also  assumed  that  no  reccmpressian  is 
necessary  after  the  diffuser  (following  heat  regeneration)  in  order  to  exhaust 
to  the  atmosphere  or  alternatively  pass  the  gases  through  a  conventional 
steam  boiler  before  rejection.  This  figure  illustrates  how  entropy  and 
ambient  enthalpy  ( temperature)  vary  through  the  cycle.  In  general  the  gas- 
velocity  will  be  low  for  the  sections  0-1,  1-2,  and  6-7-0.  Far  the  combustion 
process  2-3  no  restriction  to  subsonic  flow  has  been  assumed,  and,  in  fact, 
the  only  restriction  for  this  stage  is  that  in  practice  the  ambient  tengpera- 
ture  (enthalpy)  may  well  be  limited  (see  section  3)  • 

Fig.  1(b),  on  the  other  hand,  shows  the  same  cyole  as  Fig.  1(a)  but  here 
the  total  enthalpy  and  total  pressure  required  in  the  system  axe  shown.  In 
this  diagram  the  isentropio  sections  assumed  above  are  represented  by  a  single 
point  in  each  aase,  and  this  representation  of  the  cycle  shows  clearly  how  the 
total  pressure  in  the  system  is  reduced  progressively,  following  the  initial 
compression  in  the  leg  0-1,  back  to  atmospheric  pressure  at  0.  Also  the 
sections  of  the  cycle  are  clearly  seen  where  the  total  enthalpy,  and  hence  the 
driving  enthalpy  for  ccnveotive  heat  transfer,  is  high. 

1.1  M.H.D.  duot  temperature  and  heat  regeneration 

In  order  to  get  the  degree  of  conductivity  in  the  gas,  even  with 
seeding  ,  which  would  be  sufficient'  for  economic  use  in  the  M.H.D.  duct,  the 
ambient  temperature  in  the  duct  wouJd  have  to  be  of  the  order  of 
2500%  *  *3,4-, 5^  Although  with  st oichioms trie  mixtures  of  fossil  fuel  and  air 
at  atmospheric  temperature  it  would  appear  to  be  possible  to  obtain  tempera¬ 
tures  of  this  order  by  cunbustior.,  the  adiabatic  temperature  of  fossil  fuel 
being  around  2300-2700%,  ir.  practice  such  temperatures  are  only  obtained  by 
having  cay  gen  enrichment  of  the  air  or  using  air  preheated  (by  regeneration) 
by  about  1000°C.  Hence  when  oxygen  enrichment  proves  uneconomic  heat 
regeneration  forma  an  essential  part  of  the  cycle. 

1.2  Thrlng' a  idea  of  improved  heat  regeneration  leading  to  hypersonic1 
combustion 

The  maximum  temperature  at  whioh  a  metallic  wall  capable  of  standing  a 
very  high  pressure  difference  can  be  operated  at  the  present  time  is  about 
1001  to  1100°C.  It  is  therefore  possible  to  regenerate  heat  through  a  wall 
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maintained  at  this  temperature.  In  view  of  this  Thring  *  has  proposed  a  heat 
exchanger  which  "would  not  only  pre-heat  the  air  to  1000°C,  but  then  would 
be  extended  to  a  series  of  convergent/divergent  nozzles  in  parallel  which,  wero 
externally  heated  by  the  combustion  products,  so  that  the  heat  could  be  passed 
through  th*.  walls  and  the  expansion  to  high  velocity  could  take  place  ^ 

essentially  isothermally  at  1000°C.  In  this  way  one  arrives  at  a  gas  having 
the  necessary  velocity  and  at  1000 °C,  but  not  yet  with  the  fuel  burnt  in  it". 

Thring  goe3  on  to  point  out  that  by  this  method  one  has  transferred 
considerably  more  heat  into  the  system  at  1000°C,  due  to  the  additional  in¬ 
take  of  kinetic  energy,  than  would  be  possible,  say,  in  the  ordinary  system 
where  the  air  velocity  at  the  entry  to  the  ccmbustion  chamber  is  low.  Thus 
he  submits  that  ccmbustion  can  then  take  place  in  the  high  velocity,  possibly 
hypersonic,  stream.  Indeed  this  system  would  give  a  net  boost  to  the  total 
enthalpy  of  the  working  gas  in  the  combustion  chamber,  with  the  only  penalty, 
a  more  severe  loss  of  total  pressure  through  supersonic/hypersonic  ccmbustion 
compared  to  subsonic  combustion  (see  section  4*1),  which  would  have  to  be  & 

allowed  for  in  the  initial  ccmpression  stage  of  the  qycle.  However,  it  is 
not  clear  how  this  convergent/divergent  -  honeycomb  heat  exchanger  is  to  * 

v/ork.  For,  although  the  wall  temperature  may  be  maintained  at  1000°C  and.  the 
cooler  stream  expanded  rapidly,  heat  transfer  from  the  wall  to  the  gas  will 
be  dominated  by  the  total  enthalpy  of  the  cooler  stream  at  each  statical  rather 
than  its  ambient  enthalpy  (temperature).  Hence  the  effectiveness  of  heat- 
exchangers  of  this  type  might  be  strictly  limited  in  practice. 


2  OPEN -CYCLE  OPERATION  WITHOUT  HEAT  REGENERATION 

Although  heat  regeneration  may  be  an  essential  part  of  the  M.H.D. 
power  generation  cycle  (see  section  1.1),  in  order  to  study  in  more  detail 
the  interdependence  of  the  ccmbustion,  nozzle,  and  M.H.D.  duct  parts  of  the 
cycle,  it  is  convenient  to  include  heat  regeneration  with  combustion  heat  as 
in  Fig.  2  rather  than  to  treat  it  as  a  separate  stage  of  the  cycle  as  in  Fig.  1. 
This  can  be  justified  if  one  considers  either  the  case  of  combustion  with 
oxygen  enrichment  of  the  air  and  no  heat  regeneration,  or  more  generally  that 
for  this  simplified  model  of  the  processes  the  heat  added  by  regeneration  and 
the  heat  released  by  ccmbustion  may  be  included  under  the  heading  of 
"combustion"  alone.  Hence  for  this  latter  case  if  we  now  define  the  heat 
released  by  ccmbustion,  q,  by: 


<1  = 


%  +  ^ 


(D 

(c.f.  Figs.1  and  2) 


than,  since  the  thermal  efficiency,  n^h»  of  the  cycle  in  Fig.1  with  heat 
regeneration  la  given  by: 


^h 

with  regen. 


(2) 


the  thermal  efficiency  of  the  cycle  in  Pig. 2  considered  from  the  view  of  no 
heat  regeneration  will  be  related  by: 


^h 

with  regen* 


>th  © 

without  regen. 
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V  (V) 

without  regen.  ' 


(3) 


The  duct  effioienoy  will  of  course  be  unaltered  by  this  device  aid  is 
given  for  both  oases  by: 

^duct  =  (h  +  w  ♦  q)  =  (ht  +  w  +  qR  +  qj  (= 
o  O 

Thus  by  this  method  we  see  that  we  may  consider  the  problem  neglecting 
heat  regeneration,  but  we  must  consider  in  its  place  more  efficient  combustion 
than  is  usually  possible,  and  remember  that  the  thermal  effioienoy  of  the 
oycle  with  beat  regeneration  will  in  praotice  be  r.^jh  improved  upon  that 
without  regeneration. 

It  should  be  noted  that  in  order  to  make  the  assumptions  above,  the 
working  fluid  in  the  oyole  has  been  idealised,  and  has  been  considered  to  be 
composed  of  combustion  products  throughout  all  stages  of  the  cycle.  Henoe  a 
constant  value  for  y»  the  ratio  of  specific  heats,  for  the  working  gases  has 
been  assumed  and  the  fuel  mass  addition  has  been  neglected.  The  canbusticn 
process  (and  heat  generation)  thus  takes  the  form  of  a  latent  heat  release  in 
the  system,  or  may  be  considered  as  heat  energy  transferred  to  the  gases 
during  combust  ion* 


Energy  from  duct 
Total  energy  at 
entry  to  duct 


(4) 
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2. 1  Typical  relationships  between  the  cycle  efficiencies 

Fran  equations  (1),  (2),  (3)  and  (4)  we  see  that  for  the  cycle  given  in 


Fig. 2  we  have; 


*th 

without  regen. 


*•( 


h  +  w  +  q 
*0 


'duc  t  l  q 


For  a  typical  range  of  conditions  generally  quoted  far  this  cycle,  it 
is  shown  in  Appendix  A  that,  as  the  initial  compression  ratio,  C,  between 
states  0  and  1  increases  from  1  to  16,  equation  (5)  varies  in  the  form: 


'th 

without  regen* 


1.28. 


duct 


0.14, 


Hence  for  reasonable  duct  efficiencies  the  thermal  efficiency  in  Fig. 2 
will  for  these  conditions  be  close  to  the  duct  efficiency,  assuming  of  course 
that  the  compression,  ratio  in  turn  is  not  excessive.  However,  high  initial 
compression  ratios  may  be  requ.'  ed  if  high  pressure  losses  occur  in  the  cycle, 
(see  section  4),  and  the  above  relationship  will  no  longer  apply  (see 
Appendix  E  and  Fig.  7). 

It  is  further  shown  in  Appendix  A  that  the  thermal  efficiency  with 
regeneration  will  be  around  i«6  to  1.4  times  higher  than  the  above  value  far 
the  same  range  of  the  initial  compression  ratio. 


3  EFFECTS  OF  AMBIENT  TEMPERATURE  LIMITATIONS 
( REGENERATION V  -  NOZZIE  LEG  OF  THE  CYCLE 


THE  COMBUSTION 


As  indicated  in  section  1.1  the  total  enthalpy  of  the  system  after 
combustion  may  be  limited  in  practice,  due  to  heat  losses  in  the  combustion 
chamber.  In  order  to  simplify  the  picture  far  the  energy-exchange  processes 
in  the  combustion  chamber,  these  losses  may  be  replaced  by  considering  instead 
inefficient  or  incomplete  combustion.  However,  even  with  practical  limitations 
on  the  total  enthalpy  which  the  system  may  attain,  it  may  still  be  necessary 
to  restriot  the  ambient  enthalpy  (temperature)  throughout  the  cycle,  and  in 
particular  in  the  combustion  chamber.  This  in  general  will  not  be  the  result 
of  restricting  convective  heat  transfer,  since,  as  mentioned  earlier,  this 


7 


depends  mare  directly  on  the  total  enthalpy  of  the  system.  Nor  is  it  expected 
to  be  necessary  to  restrict  ambient  enthalpy  on  account  of  radiation  losses, 
although  radiative  energy  transfer  increases  rapidly  with  increase  of  ambient 
temperature,  since  the  ambient  temperature  ranges  considered  at  present  for 
practical  purposes  are  still  below  the  threshold  for  significant  radiative 
effects  to  occur  in  the  gas  stream.  But  a  restriction  on  ambient  temperature 
may  result  from  attempting  to  limit  or  control  the  chemistry  in  the  combustion 
reactions.  Typically  this  might  cane  fran  the  necessity  to  avoid  dissociations 
of  the  working  gas  (ionization  caning  fran  seeding)  or  influence  the  relative 
concentrations  of  caubustion  products.  In  this  latter  category,  for  example, 

L 

we  have  the  suggestion  for  inc ample te  combustion  of  fuel  in  order  to  reduce 
the  amount  of  CO ^  as  compared  to  CO,  since  at  high  temperatures  00  2  would  have 
more  serious  erosive  chemical  effect  on  the  electrodes  of  the  duct. 

3.1  Supersonic  combustion  due  to  ambient  enthalpy  limitations 

In  a  gas  stream  the  ratio  of  total  enthalpy  to  ambient  enthalpy  is 
given  approximately  by: 


i  +  be.  r  1).  n2 


Hence  the  gas  stream  is  supersonic  when: 


Uljt.il  . 

2 


(6) 


For  values  of  y  in  the  range  1.2  to  1*3  this  latter  formula  shows 
that  a  gas  stream  will  be  supersonic  when  the  total  enthalpy  of  the  stream  is 
more  than  1.1  to  1.15  times  greater  than  the  ambient  enthalpy.  Therefore,  if 
a  supersonic  M.H.D.  duct  is  to  be  used  in  the  generation  cycle  and  the  anbient 
enthalpy  level  is  restricted  during  oombustion,  these  values  for  enthalpy  ratio 
may  well  be  exceeded  at  sane  point  in  "the  combustion  chamber.  Hence  the 
combustion  prooess  may  well  be  required  to  consist  of  part  subsonic  and  part 
supersonic  combustion. 

In  Fig.  2,  since  suffix  4  denotes  entry  conditions  to  the  M.H.D. 

generating  duot,  the  inequality  (6)  above  shows  that  supersonic  combustion 

will  occur  in  the  combustion  chamber  whenever  the  maximum  ambient  enthalpy, 

h  ,  satisfies  the  inequality: 
max 


8 


^max  < 


r- 


2  +  (Y  -  1)  M 


Y  +  1 


As  mentioned  earlier  in  section  1.1,  should  correspond  to  at  least  a 
temperature  of  around  2500^.  Hence  for  y  =  1.2,  (7)  indicates  that  for 
supersonic  flow  during  combustion  would  have  oo  be  less  than  2500  ^  ■" 

(i.e.  h^)  far  =  1  and  31 80°^  (i.e.  14/1 1  h^)  for  »  2.  The  corresponding 
values  for  y  =  1*3  are  2500  \  for  =  1  and  3480 for  a*  2*  Although 
these  calculations  assume  constant  specific  heats  and  thermally  perfeot  gases, 
the  rough  values  quoted  serve  to  illustrate  hcwf  for  a  superscnio  duot,  super¬ 
sonic  combustion  may  still  be  necessary  at  some  stage,  even  when  the  ambient 
temperature  in  the  combustion  chamber  is  allowed  to  be  as  high  as  3000QK. 


4  PRESSURE  LOSSES  IN  THE  COMBUSTION  (REGENERATION)  -  NOZZTfl  LBG  OF  THE 
CYCLE 


Using  the  formula  derived  in  Appendix  B,  we  find  quite  generally  that 
the  change  (loss)  in  total  pressure  develO]«d  from  state  1  to  state  4  in  the 
cycle  illustrated  in  Pig.  2  is  given  from  (16)  by: 


Hence  for  a  given  amount  of  heat  addition  to  the  flow  we  see  that  the 
loss  in  total  pressure  has  a  lower  limit  which  depends  on  the  maximum  ambient 
enthalpy  (temperature)  that  can  be  allowed  for  the  combustion  process  (i.e. 
the  process  which  develops  the  minimum  gain  in  entropy). 

Optimisation  of  the  power  generation  cycle  therefore  requires  the 
efficiency  of  the  combustion  process  to  be  as  close  as  possible  to  this  lower 
limit  in  order  to  minimise  this  loss  in  total  pressure,  which,  as  is  seen  in 
section  2,  has  to  be  compensated  for  in  the  initial  compression  stage  0-1  of 
the  cycle. 

4* 1  Model  of  the  combustion  (regeneration)  -  nozzle  processes 

Pig. 3  shows  a  model  of  the  processes  between  states  1  (or  2)  and  4  in 
Pigs.1  and  2,  where  intermediate  states  A,  B,  C  and  D  have  been  added. 
Expansion  of  the  flow  isentropically  prior  to  combustion  (or  heat  addition)  as 
well  as  expansion  before  entry  to  the  M.H.D*.  duct  have  been  inoluded  in  this 
model. 
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Pig.  3(a)  shows  that  onoe  the  maximum  allowable  ambient  enthalpy  has 
been  reached  at  B  the  combustion  process  oould  continue  in  general  for  this 
model  to  C  at  the  same  enthalpy,  reducing  in  ambient  enthalpy  then  toD1 
before  isentropioally  expanding  to  4‘  *  However,  since  this  prooess  will  by 
the  formula  above  be  less  efficient  than  continuing  the  path  from  C  to  D  at 
maximum  ambient  enthalpy  before  isentrppio  expansion  to  4,  the  paths  to  the 
dashed  points  will  not  be  considered,  and  the  point  C  may  be  Ignored  in  the 
model  of  the  cycle  shown. 

Appendix  C  gives  the  relationships  required  to  determine  the  various 
states  through  the  isentropic  nozzle  prooesses,  1  (or  2)  -  A  and  D  -  4,  and 
the  combustion  prooesses  A-B  and  B-D.  In  the  latter  category  we  have  constant 
antoient  enthalpy  o embus ti on  for  B-D.  The  combustion  prooess  A-B  is,  however, 
not  defined  and  for  this  one  may  assume,  as  in  Appendix  C,  constant  pressure 
(velocity),  constant  Mach  number,  constant  area  (which  would  Include 
combustion  through  a  detonation  wave) ,  or  seme  mare  general  combustion 
process,  determining  properties  in  this  latter  case  by  using  the  basic 
formulae  of  Appendix  B. 

Comparative  study  of  the  loss  of  total  pressure  required  by  the 
combustion-plus-nozzles  model  illustrated  in  Pig.  3  can  thus  be  made  by 
defining  (which  in  sane  cases  may  be  taken  to  be  so  large  that  in  fact 

no  ambient  temperature  restriction  exists)  and  varying  both  the  combustion 
process  A-B  and  the  velocity  at  state  A.  The  velocity  at  state  4,  the 
conditions  at  state  1  (or  2),  and  the  total  combustion  heat  release  being 
fixed  for  this  investigation,  (the  parameter  A  being  fixed  automatically 
by  the  other  assumptions) .  Again  from  the  formula  for  the  overall  loss  of 
total  pressure  given  earlier  in  this  section,  it  is  obvious  that  the  minimum 
loss  of  total  pressure  (minimum  gain  in  entropy)  will  result  whan  the  combus¬ 
tion  path  A-B  is  chosen  such  that  the  ambient  enthalpy  (temperature)  at  each 
stage  of  the  process  is  as  large  as  possible  (i.e.  as  close  to  hmnY  as 
possible)  . 

4«  2  Simplified  calculations  based  on  the  above  model  allowing  for  completely 

supersonic  and  completely  subsonic  combustion 

Pig. 4  shows  the  results  from  calculations  (made  by  Thompson)  based  on 
the  model  in  Fig. 3  for  a  fixed  amount  of  heat  addition  by  combustion  where 
the  combustion-nozzle  processes  between  states  2  and  4  have  been  considered, 
(i.e.  only  the  processes  upstream  of  the  M.H.D.  duct  which  occur  after  the 
heat  regeneration  leg  1-2  have  been  considered.)  AI30  no  restriction  on 
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ambient  enthalpy  level  In  the  combustion  process  has  been  made*  Both 
completely  eijperseniQ  and  completely  subsaiio  combustion  processes  have  been 
considered  together  with  the  appropriate  iaen tropic  nozzle  expansions  which 
are  required  to  match  in  with  these  in  order  to  arrive  at  the  same  Mach 
number,  M^  *  1*4,  at  entry  to  the  M.H.D.  duct. 

The  key  to  the  various  oombustion-nozzle  combinations  studied  is  given 
in  Fig*4(ft),  with  the  various  paths  2-I^-4  identified  by  the  suffix,  (n), 
of  the  intermediate  state  1^ .  For  the  ocnditicns  of  this  example  it  was 
not  possible  to  include  the  path  with  isentrcpic  expansion  followed  by  super¬ 
sonic  combustion  at  constant  area.  The  reason  far  this  is  that  in  this  case 
the  Initial  total  enthalpy  at  state  2  was  not  sufficiently  high  to  allow 
expansion  of  the  flow  to  the  large  velocity  required  for  matohing  this  process 
to  constant  area  supersonic  o embus, tion  through  to  a  Mach  number  of  1*4  at 
state  4* 

Fig. 4(b)  illustrates  clearly  for  this  example  the  large  loss  of  total 
pressure  in  the  system  through  such  supersonic  combustion  processes.  As 
mentioned  earlier,  such  pressure  losses  would  thus  have  to  be  tolerated  in  the 
Thring  model  of  the  heat  regeneration  -  combustion  system  described  in 
section  1.2. 

Fig. 5  shows  results  from  calculations  similar  to  those  illustrated  in 
Fig. 4,  where  the  Maah  number  at  entry  to  the  duct  was  assumed  to  be  2  instead 
of  1.4*  Again  matohing  supersonic  combustion  at  constant  area  to  the  required 
isentropic  nozzle  expansion  was  not  possible  for  the  conditions  assumed.  Nor 
was  this  matching  ix,  fact  possible  for  supers cnio  combustion  at  constant 
pressure,  although  apart  frem  this  the  general  picture  is  much  the  same  as 
before  in  Fig.4» 

The  processes  shown  in  Figs*.'r  and  5  are  by  no  means  optimum,  but  they 
illustrate  the  large  pressure  losses  involved  with  having  fully  supersonic 
combustion  in  the  cycle  as  compared  to  fully  subsonic  combustion. 

5  CONSTANT  MACH  NUMBER  M.H.D.  DUCT  -  INTEGRATION  OF  THE  DUCT  INTO  THE  CYHT.tc 

In  order  to  estimate  the  efficiency  of  the  M.H.D.  duct  itself,  the  duct 
with  constant  Mach  number  has  been  considered.  In  fact,  as  is  shown  in 
Refs. 6  and  7,  this  daot  is  close  to  ha/ing  optimum  efficiency,  even  when  a 
constant  load  factor,  K  (see  Appendix  D) ,  is  assumed. 
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A  summary  of  the  relations  appropriate  to  this  duct  are  given  in 
Appendix  D,  and  the  efficiency,  ri^  j_,  far  such  a  system  is  illustrated  in 

U.UC  b 

Pig. 6.  In  this  figure  two  duct  Mach  numbers  have  been  chosen  together  with 
two  values  of  the  load  factor,  K.  It  should  be  noted  hew  sharply  (almost 
exponentially)  the  pressure  ratio  required  aoross  the  duct  rises  with  a 
demand  far  increase  in  efficiency.  In  order  to  integrate  the  duct  into  a 
complete  cycle  it  is  therefore  obvious  that  a  large  total  pressure  loss  will 
have  to  be  allowed  for  in  the  duct  itself,  and  this  means  that  loss  of  total 
pressure  will  thus  be  at  a  high  premium  throughout  the  early  part  of  the  cycle 
upstream  of  the  duct. 

5*1  Integration  of  the  constant  Mach  number  duct  into  the  complete  cycle 

As  is  shown  in  Appendix  D  equation  (24)  the  M.H.D.  duct  efficiency- 
possible  in  any  particular  case  may  be  fixed  by  defining  the  ratio  of  total 
pressures  (or  static  pressures)  across  the  duct.  In  practice  the  duct  length 
required  to  achieve  this  efficiency  might  be  too  long  if  the  temperature  level 
or  magnetic  flux  strength  in  the  duct  is  limited.  However,  in  order  to  obtain 
an  idea  of  the  relationship  of  overall  cycle  efficiency  to  duct  efficiency  it 
will  be  assumed  that  the  duct  efficiency  defined  by  the  ratio  of  total 
pressures  across  the  duct  is  realistic. 

When  this  model  of  an  M.H.D.  duct  system  is  integrated  into  the  full 
thermal  cycle  the  overall  thermal  efficiency  is  given  by  equation  (28)  in 
Appendix  E  asi 


th 

without  regen. 


'duct 


1  +  £ 
)  Y 


-PC, 


This  equation  illustrates  that  fer  a  fixed  amount  of  heat  addition  by 
combustion,  q,  the  overall  thermal  efficiency  decreases  as  entropy  increases 
due  to  less  favourable  choice  for  the  combustion  process. 

When  the  type  of  combustion  process  is  fixed,  the  relationship  of  overall 
thermal  efficiency  to  M.H.D.  duct  efficiency  for  various  total  amounts  of  heat 
addition  during  combustion  has  also  been  studied  in  Appendix  E.  Explicit 
optimisation  was  not  found  to  be  possible  far  most  of  the  combustion  processes 
considered,  but  for  the  case  of  constant  pressure  (velocity)  combustion  at  low 
subsonic  speeds  it  was  found  that  the  optimum  value  for  the  thermal  efficiency 


.,,**«*>  WWi>Wr,J3  J^WPU.tTj  j/maMfefttert 


12 


is  obtained  by  having  the  largest  amount  of  heat  addition  possible  during  the 
combustion  process*  Further  the  same  trend  is  true  for  this  oase  at  higher 

y2 

flow  velocities  ^given  by  the  velocity  parameter  h^)^ 

o 

chamber,  as  can  be  seen  from  the  examples  illustrated  in  Pi9*7{e)»  (b)  and  (o). 
Here  the  cycle  efficiency  is  plotted  against  the  combustion  total-heat-addition 
parameter  and  is  shown  by  solid  curves  for  each  value  of  the  velocity 

parameter*  The  regions  in  these  figures  for  fully  subsonic,  fully  supersonic, 
and  part  subsonic  and  part  super  sen  ic  combustion  are  illustrated  by  the  dashed 
lines.  (The  Mach  numbers  assumed  for  the  M.H.D.  duct  in  these  examples  were 
1*4,  0.5  and  2*)  Thus  it  is  seen  that  the  M.H.D*.  duct  efficiency  of  0*2 5  is 
approached  only  for  the  case  of  subsonic  combustion  with  high  heat  addition, 
the  supersonic  combustion  region  being  much  less  efficient  due  to  the  large 
initial  compression  which  is  required  to  canplete  the  oyole  for  this  type  of 
combustion  process* 


6  CONCLUDING  BMf  ABKH 

In  the  above  simplified  examination  of  parts  of  the  M.H.D*  power 
generation  cycle  no  allowances  for  heat  loss  or  real  gas  effects  have  been 
made  in  the  analysis.  Although  it  might  be  argued  that  efficient  regeneration 
processes  should  compensate  to  some  extent  far  heat  losses,  the  design  of  a 
system  which  requires  the  minimum  amount  of  regeneration  in  this  way  might 
well  prove  to  be  necessary  in  practice. 


Thus  far  the  operation  of  an  M.H.D.  paver  generation  duct  of  supersonic 
speeds,  the  possibility  of  having  an  effectively  ’cold*  nozsle,  perhaps  with 
regeneration  as  Thring  suggests,  followed  by  supersonic  combustion  might  be 
desirable,  since  then  no  hot  nozzle  expansion  is  necessary,  the  ambient 
temperature  is  limited  (see  section  3),  and  the  maximum  total  enthalpy  and 
ambient  temperature  are  produced  only  at  the  entry  to  the  M.H.D.  duct. 
However,  as  is  shown  in  section  4*2,  pressure  losses  may  prove  to  be  too  hi^ti 
in  such  a  system,  and  because  of  the  large  initial  compression  required  to 
compensate  for  this  the  overall  efficiency  may  be  low,  as  was  indicated  from 
the  examples  studied  far  supersonic  combustion  in  section  5* 


Since,  on  the  other  hand,  canplete  combustion  before  expansion  in  a 
nozzle  to  the  required  conditions  for  the  M.H.D.  duct  produces  problems  of 
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Appendix  A 

TYPICAL  RELATIONSHIPS  BETWEEN  THE  CYCLE  3ZPICIENCIES 


In  section  2.1  we  have: 
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For  typical  cycles  shown  in  Figs.1  and  2  a  rough  estimate  of  the  order 
of  magnitude  generally  quoted  for  each  of  the  quantities  shown  is: 

ht  -  300°K 
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=  5  4 


(-  1500°K) 
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w  =  (C) Y  -  1  h 
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where  C  is  the  initial  compression  ratio  fran  state  0  to  state  1, 
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(based  on  the  fact  that  the  limit  of  regeneration  is  given  by: 

ht  =  ht  +  t  +  =  4  ht  1 200^) . 
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Hence: 
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For  a  typical  leg  of  die  qycle  «-83  where  an  amount  of  heat,  q,  is 
either  added  to  or  produced  by  combustion  in  the  flow,  we  have  the  following 
relations  when  an  ideal  gas-model  is  assumed; 


From  (9),  (10)  and  (ll)  we  have: 


Energy  equation  (8) 


Moment  vim  equation  (9) 
Equation  of  state  (10) 


Definition  of  h^  (11) 


Definition  of  (12) 
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Using  (10)  and  (20)  we  obtain  a  quadratic  for  vi 


(Y  +  1)  v  -  2r  v 


f  — +  A»  +  2  (y  -  1)  (ht  +  q)  ■  o 
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The  realistic  solution  for  v  fran  this  quadratio  is: 
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...  (21) 


Thus  using  (21)  together  with  (20),  (11),  (12)  end  (17)  the  properties 
at  point  0  may  be  determined.  It  should  be  noted  that  in  this  oase  it  is  not 


convenient  to  calculate 


•—L  in  order  to  determine  p  directly.  Note 
h  *0 


from  (8),  (11),  (12)  and  (20)  we  have: 
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FOIthiULAE  RELATING  THE  VARIOUS  PARTS  OF  THE  MODEL  FOR  THE 
OOMBU  ST  ION  -PT  JiR  -N0Z7.T  ftfl  SECTION  OF  THE  CYCLE 
(see  Fig* 3) 

(i)  The  process  fran  1  (or  2)  to  A  in  Fig* 3  is  an  isentropic  nozzle 
expansion  and  is  given  by: 


ht  "  ht 
A  1 


(°r  h  ) 
t2 


\  =  % 


(or  P+  ) 
*2 


v  (<  (2  h  )z)  depends  upon  the  nozzle  expansion,  and  oan  therefore  be  chosen 
A  t* 


arbitrarily. 


(ii)  Process  from  A  to  B  is  given  by;  (h^  s  h 


h,  =  h  +  Aq 


together  with  (for  the  following  special  cases  given  in  Appendix  B)  i 
(a)  Constant  static  pressure  (velocity) 
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(b)  Constant  Mach  number 
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and 


p  2=± 

Ay 


(e)  Constant  area 
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^iii)  Process  fran  B  to  D  is  given  by:  (h^  s  h 

h  =  h  +  (1  -  X)q 

VD  ■  fvv 

(or  tip  =  hs) 

and  from  the  constant  ambient  enthalpy  process  given  in  Appendix  B 
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(iv)  Process  from  D  to  4  is  again  an  isentropic  nozzle  expansion  given  by: 


P*  = 


with  v,  (<  V2  h7)  chosen  to  give  required  velocity  at  state  4* 
4  t4 
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FORMULAE  RELATING-  HtOESRTIBS  IN  A  CONSTANT 
MACH  NUMBER  U.H.D.  DUCT  . 

Fran  Ref. 6  using  the  notation  of  Fig.  2  we  find  that  for  the  constant 
Mach  number  M.H.D.  duct  the  efficiency  of  the  duct,  ‘nduct*  and  temperature 
and  pressure  ratios  across  it  are  given  by: 


where  P  =  [l  +  ^  (1-K)  (y-l)  M^J  y/(y-1)  K,  and  M^  and  y  are  constant. 

K  is  considered  oonstant  in  this  analysis  (0  <  K  <  1)  and  is  the  external 

load  factor  .(■  h  <see  Ust  oe  ^°^>)  «*  duct-  A  ‘^.1 

value  fa:  K  being  0*5  to  0.75  far  an  efficient  duct. 

Fran  (22)  and  (23)  we  obtain: 

(’-W0  • 

This  relationship  is  illustrated  in  Fig.  6  arxl  shews  how  the  ratio  of 

total  pressure  across  the  M.H.D.  duct  varies  with  the  efficiency  of  the  duct, 

for  v  =  1*3  and  various  values  of  M.  and  K. 

4 

Although  for  particular  values  of  M^,  y  and  K  the  duct  efficiency 
depends  only  on  pressure  and  temperature  ratios  across  the  duct,  the  duct 
length  required  to  achieve  this  efficiency  depends  on  the  absolute  values  of 
the  conditions  in  the  duct.  It  must  always  be  remembered,  therefore,  when  a 
value  for  the  duct  efficiency  is  assumed,  that  in  practice  this  value  may 
correspond  to  an  unrealistic  length  for  the  duct,  when,  for  example,  the 
temperature  level  or  magnetic  flux  strength  is  limited.  Also  it  should  be 
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noted  the!  optimum  duct  efficiency  does  not  necessarily  correspond  to  optimum 
overall  efficiency  for  the  full  thermal  oycle. 

For  a  constant  Mach  number  duct  with  constant  K  and  y  the  duot  length 
x  is  given  by* 


x 

X* 


1  -  (1 


T'duct> 


where  w'  =  (3  -  ■?)  -  (y*  -  Pa1)# 


x* 


V  \ 

K(1  -  K)  V  .  <r4 


and  the  indices  y*  and  z'  originate  from  an  expression  of  the  form  O’  *  lr  p 
for  the  electrical  conductivity  cr.  Typically  y*  is  of  the  order  of  IQ-13  end 
z'  is  of  the  order  of  0«4-0»5* 
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INTEGRATION  OF  THE  CONSTANT  LLACH  NUMBER 
M.H.D.  DUCT  INTO  THE  COMPLETE  CYCLE 


Using  the  notation  in  Pigs. 2  and  3*  Tran  Appendix  B  equation  (16)  we 
have  direotly,  by  considering  the  energy  addition  w  during  isentropic 
compression  and  heat  addition  q  during  canbustion: 


q  +  w  + 


q  +  w  + 


Jo  h 


where  As  is  the  entropy  gain  in  the  process. 

If  we  assume  that  p  =  p.  (i.e.  that  the  duct  operates  down  to  the 

initial  pressure) ,  then  by  matching  the  pressure  through  the  system  we  see 
that: 


Prom  equation  (26)  and  (24)  in  Appendix  D  it  follows  that: 


This  shows  how  an  increase  in  initial  compression,  given  by  the  energy 
parameter  w,  arises  from  an  increase  in  entropy  in  the  system. 
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Prom  (5)  Beotian  2.1,  (21*)  Appendix  D,  and  (27)  eliminating  w  we  get: 


^th 

without  regen* 
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[(1  -  r'iuct) 


Equation  (28)  illustrates  how  far  a  fixed  amount  of  heat  addition,  q, 
the  themal  efficiency  of  the  qyole  decreases  as  entropy  in  the  system 
Increases*  Hence  the  most  efficient  system  for  a  fixed  value  of  q  requires 
the  minimum  entropy  gain  or  total  pressure  loss  in  the  oombustion  part  of  the 
cycle. 

This  is  an  obvious  result  and  is  associated  with  the  comparison  of 
various  combustion  processes  for  a  fixed  amount  of  heat  addition,  q*  However, 
for  a  particular  combustion  process,  optimisation  of  the  thermal  efficiency 
when  q  is  varied  is  far  more  complex,  since  in  (28)  As  is  a  function  of  q. . 

Differentiating  (28)  with  respect  to  q  it  follows  that  for  stationary 


values  of  tj. 


,  q  must  satisfy  the  equation: 


without  regen. 


fcfr]  - 


~  1  "r) 


w* 


Hence  if  a  clear  maximum  value  of  exists  it  will 

without  regen* 

correspond  to  a  value  far  q  given  by  (29).  In  (29)  As  may  be  assumed  to  be 
a  function  of  q  and  is  found  from  (27),  assuming  w  also  as  a -function  of 
q,  and  the  appropriate  law  for  the  combustion  process  assumed* 
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Appendix  £ 


As  examples  of  this  optimisation  process  we  may  consider  in  detail  the 
cases  of  ( i)  constant  static  enthalpy,  (ii)  constant  static  pressure 
(velocity),  and  (iii)  constant  Viach  number,  combustion,  which  are  listed  in 
A  pendix  B.  For  these  processes  we  may  consider  combustion,  using  the 
notation  of  Figs.l  and  2,  to  start  from  some  state  A  between  1  (or  2)  and  4 
as  in  Fig.  3» 

Hence  we  have  in  addition  to  equation  (27)  fcr  each  case  the  following 
equations  from  Appendix  Bs 

m  IS  =  -3- 


(w  ■f  h.  -  £  v*) 
o  A 


(q  +  h  -  i  vf) 

A  A 

K  "  *  'I' 


(q  +  w  ♦  ht  -  £  vp 
o 

(  ’"  ♦  ht  -  i  vp 

O 


q  +  h 


(iii) 


_  (  _ k\ 

-  \  V  ) 


!i'\  ehA 


q  *  w  ♦  h. 


w  +  h 


w  ♦  hL  -  2  v 
>o  A 


V  w  +  ht  ) 


Appendix  B 


Particular  gage: 


The  simplest  case  above  to  study  algebraically  is  that  of  constant 
static  pressure  combustion*  For  this  case  equations  (27)  and  (30)  (11)  ®ay 

is 

C  ' 

be  solved  f or  e  p  and  er  direotly  to  give: 


‘(t) 

o 

-  ^(r-1)  /  \  \  /  q  \  "T 

-  w  -Wr) 


-  %r-i)  /  \  /  a  \ 

[<’  -  Vct>  Y  (32) 

*  A  A 


whore 


9  ■  [[•h'*(7^*)+  (1  't,4>“^y  T*4(r^r)<1  ”ltao^ 


-  -®(y“1  ! 

’i  Y 


W 

However,  substitution  of  e  p  from  equation  (3^)  into  (29)  results  in  an 
equation  for  q  which  is  not  amenable  to  simple  analytical  treatment*  Hence 
it  is  better,  to  stuiy  the  behaviour  of  the  thermal  efficiency  with  q  more 


directly  by  substituting  for  e  p  from  (3"0  into  equation  (28).  This  gives: 


'th 

without  regen. 


1  J1'*  y  "  W 

■  1  ‘  TT^V1  ',1auot)  _  1 


/  \  1  +  ^  -  0 

y 

-  ^(y-1)  /  \  /  , 


w Y 


-  -  (-A  ♦  © 

V2hJ  \htJ  J 


...  (33) 
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Apptadix  E 


At  any  stage  during  the  combustion  process  we  have: 


,•29 


v?  *  (y  -  1)  MZ  h  s  (y  -  1)  M*  (w  +  h.  ♦  q'  -  4  v*) 


(5k) 


where  q'  ie  the  amount  of  heat  added  up  to  that  stage  during  ecnbustion  and 
0  <  q*  <  q» 

Hence ,  frem  (34),  the  flow  during  ecnbustion  will  be  sonio  cr  super- 
sonio  at  points  where * 


TA 


>  (Y  -  1)  (w  +  ht  +  q'  -  ?  V*) 


i.  e. 


i  ffi)  \  *  \  ♦»*>  • 


[v>) 


Thus  since  q’  lies  in  the  range  0  <  q*  <  q,  the  ecnbustion  process  vrfil  be 
completely  supersonic  when: 


1 

2  v  - ') 


+  1  \  2  /  «  \ 
vA  >  (w  +  ht  +  q) 

o 


and  canpletely  subsonic  when: 


1  /XJ lX\ 

2  Vy-  V 


vA  <  (w  +  ht  ) 
o 


U3 


(36)  and  (32)  may  therefore  be  used  to  define  the  regions  for  completely 
supersonic  and  completely  subsonic  combustion.  In  between  we  have  a  region 
where,  in  this  oonstant  velocity,  v^,  case,  the  flow  starts  off  supersonic  but 
gradually  reduces  to  subsonic  as  heat  is  added  during  combustion. 

As  a  special  case  of  equation  (33)  we  have  the  case  of  low  subsonic 


0r) 


combustion  when  [-■"?—  J  may  be  negleoted.  Hence  we  get  in  this  oase: 

'o 


r,th 

without  regen. 


=  n 


duct 


-  v=t>  Y  -1]  • 


S3 


1 


1 
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Appendix  £ 

Equation  (37)  shows  that  for  oonbustion  at  low  subsonic  speeds  tbs  ' 
maximum  value  of  the  thermal  efficiency  is  obtained  for  the  largest  value  of 

^  *  ^Slnue  >  1,  see  (23)  Appendix  D.J 


In  the  oa;,v  of  higher  flew  velocities  in  the  combustion  chamber, 
equation  33)  in  faot  gives  a  similar  trend  for  the  variation  of  the  thermal 
efficiency  with  the  total  amount  of  heat  addition  during  combustion*  This  is 
cJ  y  seen  in  Fig* 7(a),  (b)  and  (c)  where  oases  associated  with  a  constant 
Vauh  uum.ar  M.H.b.  duct  (K  =  0*75) »  for  Mach  numbers  of  1*4,  0*5  and  2 
respectively  and  duct  efficiency  of  0.23,  a re  plotted  far  various  values  of 


the  oonbustion  velocity  parameter 


•  Also  shown  are  the  regions  far 


fully  subsonic  and  fully  supersonic  flow  during  combustion  which  are  given  by 
equations  (36)  and  (32). 


31 


f 


f 

i 


*2 


saiBOia 

T  temperature 

v  velocity 

h  enthalpy 

t  jtal  enthalpy 

M  Mach  number 

heat  added  by  regeneration 

q^,  heat  added  by  combustion 

q  total  heat  addition 

q*  intermediate  heat  addition 


P 

pt 

P 

Y 

C 

e 

A 

B 

X 

x* 


y' 

z' 

w' 

3 

or 

K 


heat  addition  on  going  from  state  a  to  state  0 
specific  heat  at  constant  pressure 
entropy 

^Defined  by  ~  =  J  — j*- 

energy  added  during  initial  isen tropic  compression 
ambient  pressure  (N/m  or 

total  pressure  (N/m  or 

energy  output  from  1.1.  H.D,  duct 
ratio  of  specific  heats 
initial  compression  ratio 

expression  defined  in  Appendix  D  for  equations  (31)  and  (32) 
expressions  defined  in  Appendix  A 
Expression  defined  in  Appendix  A 
length  of  M.H.D.  duct 

'•characteristic"  length  far  M.H.D,  duct  given  by 
equation  (25)  in  Appendix  D 

index  used  in  Appendix  D  for  or 

index  used  in  Appendix  D  for  <r 

index  used  in  Appendix  D  equation  (25) 

index  used  in  Appendix  D  equation  (23) 

electrical  conductivity 

external  load  factor  for  M.H.D.  duct 


(n/soc) 

(m2/eeo2) 

(m2/aec2) 

(m2/sec2) 

(m2/seo2) 

(m2/seo2) 

(m2/seo2) 

(m2/seo2) 

(m2/sec2°K) 

(mt'/sect'k) 

(m2/  sec2) 
atmospheres) 
atmospheres) 

(m2/sec2) 


(m) 

(m) 
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SYMBOLS  (Con tdO 

S  electric  field  across  the  duct 

B  magnetic  field  across  the  duct 

Suffices 


t  total  conditions 

stages  in  the  cycle 


(webera/  m  sec  a  volta/m) 

2 

(webera/m  ) 
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00  SUBSONIC  COMBUSTION  AT  CONSTANT 
AREA  —  ISENTROPIC  EXPANSION 
(n)  J  (iii)  ISENTROPIC  EXPANSION -SUPERSONC 
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COMBUSTION 


NOTATION  FOR  COMBUSTION 
PROCESSES  2-1  (h)-  4, 


(n)  AS  ABOVE 
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b  TOTAL  ENTHALPY  v  TOTAL  PRESSURE  (atmos) 

FIG.5a*b  RESULTS  FROM  MODEL  GIVEN  IN  FIG.  3  WITH 

COMPLETELY  SUPERSONIC  AND  COMPLETELY  SUBSONIC 

COMBUSTION  PROCESSES  .  MACH  NUMBER  AT  ENTRY  TO 

M.H.D.  DUCT  M4=2,  ff  =  l3  (AFTER  THOMPSON) 
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FIG.  6  EFFICIENCY  OF  A  CONSTANT  MACH  NUMBER 
M.H.D.  DUCT  VERSUS  RATIO  OF  TOTAL  PRESSURE 
ACROSS  THE  DUCT  (8  =  I  3) 


CONSTANT  PRESSURE  (VELOCITY)  COMBUSTION  AND  A 
CONSTANT  MACH  NUMBER  M.H.D.  DUCT  (K=075)  OF 


EFFICIENCY  0-25,  FOR  VARIOUS  VALUES  OF  THE 


COMBUSTION  VELOCITY  PARAMETER 

(see  APPENDIX  E  (33)  -  (37)) 


va‘ 
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